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Abstract We have employed the improved Peierls-Nabarro (P-N) equation to study the properties of 1/2(1 10) edge dislo- 
cation in {1 1 1} plane in FCC metals Al, Cu, Ir, Pd, and Pt. The generalized-stacking-fault energy (GSFE) surface entering 
the equation is calculated by using first-principles density functional theory (DFT). The accuracy of the method has been 
tested by calculating values for various stacking fault energies which favorably compare with the previous theoretical and 
experimental results. The core structures, including the core widths both of the edge and screw components, dissociation 
behavior for edge dislocations have been investigated. The dissociated distance between two partials for Al in our calcu- 
lation agrees well with the values obtained from the numeric simulation with DFT and molecular dynamics simulation, as 
well as experiment. Our calculations show that it is preferred to create partial dislocation in Cu, and to be easily observed 
full dislocation in Al, Ir, Pd, and especially Pt. 
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1 Introduction 

It is widely accepted that the dislocations have a direct influence on the deformation mechanism of materials 

IQ , and there has been a great deal of interest in describing accurately the dislocation core structure on an atomic 

scale because of its important role in many phenomena of crystal plasticity (2l|3l. In FCC metals dislocations 

can reduce elastic energy by dissociating into Shockley partial dislocations connected by a stacking fault HI. 

But as the core regions of two partials overlap, whether the partials can be observed has not been governed by 

the stacking fault energy but also generally by atomic interactions. The dislocations in Al have always been 

focused on in the past twenty years @[5j[6ll3|8j|9j[T0im[T2l[T2l]. The accurate prediction of the atomic-scale 

core geometries for dislocations in Al is needed to inform higher-length scale models of plasticity in Al alloys, 
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as well as as other FCC metals [4]. In addition, it is well known that the nature of slip cannot be described 
in terms of an absolute value of stacking fault energy and a correct interpretation requires the generalized- 
stacking-fault energy (GSFE) curve, involving both intrinsic (stable) and unstable stacking fault energies [ 14]. 

There are usually two types of theoretical approaches which have been employed to study the core proper- 
ties of dislocations. The first method is based on direct atomistic simulations or first-principles calculations, and 
the second method is based on the the framework (Peierls-Nabarro) P-N model. Empirical interatomic poten- 
tials involve the fitting of parameters to a predetermined database and hence may not reliable in describing the 
core properties. On the other hand, first-principles calculations, though considerably more accurate, are compu- 
tationally expensive for studies of dislocation properties. The P-N model, in which the discreteness of crystal is 
neglected, is essentially a continuum treatment, but the dislocation core, the region of inelastic displacement, is 
given an approximate atomistic description. In addition, the partials in FCC metals are mixed dislocations with 
both edge and screw components H. Schoeck and Mryasov et al adopted the generalized 2D P-N model to 
consideration the two-component displacement field of dissociated dislocation |[8][T5][T6]|, but discreteness of 
crystals is neglected. Recently, the 2D improved P-N dislocation equation taken into account the discreteness 
of crystals based on the lattice dynamics and the symmetry principle has been obtained lfT71[T8"ll . Both for the 
P-N model and improved P-N equation [18], the forces in the dislocation core, where the atomic atomic-scale 
discreteness really counts, are currently approximated with the GSFE surface (y-surface) ||3j[T3j[l9l|20l. The 
GSFE surface is the interplanar potential energy, which can currently be calculated accurately from ab initio 
calculations, for sliding one half of a crystal over the other half. 

In this paper, we carry out first principles calculations on the GSFE surfaces (y-surface) for FCC metals 
Cu, Al, Ir, Pd, and Pt. The y-surfaces are applied in the improved P-N equation to the edge dislocations. The 
accuracy of the method has been tested by calculating values for various stacking fault energies which favorably 
compare with the previous theoretical and experimental results. Dislocation structures are calculated, including 
the core widthes and dissociated distance between two partials. This paper is organized as follows. In section 
2, we carry out the first principles calculation of GSFE surface on (111) plane in FCC metals Cu, Al, Ir, Pd, 
and Pt. In section 3, the core structure of dissociated dislocation is obtained from the improved P-N equation 
with variational method, and dislocation properties have also been discussed. 

2 First principles calculations of generalized-stacking-fault-energy surface 

2.1. Computational Methodology 

In the present study total-energy calculations based on the density functional theory (DFT) embodied in 
the Vienna ab simulation package (VASP) ETl 1221 1231 is employed. The Perdew-Burke-Ernzerhof (PBE) 



2 



E4l l25ll exchange-correlation functional for the generalized-gradient-approximation(GGA) is used. A plane- 
wave basis set is employed within the framework of the projector augmented wave (PAW) method [26 ,(271. On 
the basis of our tests, we have chosen the plane-wave energy cutoff of 500 eV for all calculated metals. An 
initial calculation is undertaken to determine the optimum lattice constants as well as the elastic constants. For 
the first-Brillouin-zone integrals, reciprocal space is represented by Monkhorst-Pack-special k-point scheme 
[28] with 15 x 15 x 15 grid meshes for these initial calculations, while the generalized-stacking-fault energy 
(GSFE) calculations employ 21 x 21 x 3 grid meshes for our calculated FCC metals Cu, Al, Ir, Pd, and Pt. 
The equilibrium theoretical lattice structures are determined by minimizing the Hellmann-Feynman force on 
the atoms and the stress on the unit cell. The convergence of energy and force are set to 1.0 x 10~ 6 eV and 
1.0 x \0~ 4 eV/A, respectively. In Table 1, we give the all equilibrium lattice constants and elastic modulus for 
all studied FCC metals in our calculation, and it is shown that the results agree well with the experiment values 
MM- 

2.1. The generalized-stacking-fault energy surfaces 

In the present study, we calculate the GSFE surfaces in the closed-packed (111) surface which is the trigonal 
lattice, since the slip between the closed-packed surfaces is most easily for FFC metals. The ideal FCC struc- 
tures have the configuration • • • ABC ABC ■ ■ • stacking sequence of the atomic planes. To simulate the block 
shearing process we use a slab consisting of twelve atomic layers in the (111) direction. Between periodically 
repeated slab the vacuum space of 15A normal to (1 1 1) plane is chosen to avoid interactions between two slabs. 
The stacking-fault energies (SFEs) for the slab consisting of twelve and fourteen atomic layers are reasonably 
close. In addition, we find that the fluctuations of calculated results for vacuum gap of 15 A and 18A are less 
than 1.0%. So, the adequate convergence with respect to the supercell size consisting of twelve atomic-layers 
thick slab and vacuum gap of 15A is indicated. The slab calculations is appropriate since the fault interactions 
are short ranged [13]. Comparison with a six atomic-layers thick slab calculations of intrinsic-stacking-fault 
(ISF) by Hartford et al |[T3l , our twelve atomic-layers thick slab calculation is more accurate, but with great 
computational consumption. 

The SFEs are determined as a difference of total energies for two sides of the slab, in which there are six 
atomic layers for below and above, designed to simulate faults with vectors u = and u t 0. The stacking- 
fault (SF) vector u, which is a two-dimension (2D) vector (u x , u y ), is the displacement of the upper half-crystal 
relative to the one below. The ISF configuration corresponds to a slip of A= in the (112) direction, resulting in 
the SF vector u = ^(-^,5) (where a is the lattice constant). The unstable stacking-fault energy (USFE) yusF 
corresponds to the lowest energy barrier that needs to be crossed for the slip from the ideal configuration to 
the ISF configuration in the (112) direction. In order to obtain the correct SFEs, structural optimizations must 
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be considered in this study. The relaxations perpendicular to the (111) plane are carried out by a combined 
optimization of supercell volume and atomic coordinates in shifted configuration. 

To obtain the GSFE surfaces of FCC metals, we have calculated SFEs for the slips along the (112) and 
(110) directions, respectively. Taking Al as an example, the fitted GSFE surface is shown in Fig. 1, while the 
GSFE surfaces for other FCC metals has the similar shapes. The first energy maximum encountered along the 
(112) direction is the USFE yusF> which represents the lowest energy barrier for dislocation nucleation[?]; the 
first energy minimum value at 1/6(112} (a/ V6) corresponds to the ISF configuration, where a full dislocation 
dissociates into a pair of Shockley partials. The projection of the GSFE surface on the (110) direction are 
symmetric with respect to the slip displacement of V2a/4. In Fig. 2, we show the projection of the GSF energy 
surface on the (112) (Fig. 3(a)) and (110) (Fig. 3(b)) directions. The values of unstable stacking-fault (SF) 
energy along the (110) direction are found to be larger than one along the (112) direction. The energy values 
of the various stacking faults obtained from the DFT calculations are summarized in Table 2. We also compare 
our results with the previous theoretical results lTT3l [3111321 and experimental results |[T][33j|34]]. The calculated 
results of ISF energies (ISFEs) yisF are in good agreement with the experimental values. Comparison with the 
results obtained from the EAM (embedded atom method) (3JJ or MEAM (modified EAM)|[32l by Zhang et al, 
our results agree better with the experimental results. Our results for Al and Pd are also in consistent with the 
results obtained from Hartford et al [Q2] by first-principles calculations. The value of yj for Pt is 282mJ/m 2 
which is in better agreement with the experimental value 322mJ/m 2 than the results lllmJ/m 2 obtained from 
MEAM 021 . 

The GSFE surface y(u) satisfies the translational symmetry, 

y(R + u) = y(u), (1) 

where R is the lattice vector. The GSFE surface can be represented with the aid of reciprocal lattice by a 2D 
Fourier series which reflects the full symmetry of the (111) plane, 

y(u) = Y J 7Ge iGu (2) 
G 

where G is the reciprocal lattice vector. While the nearest neighbor approximation in reciprocal lattice space is 
taken into account 

i /4nu y \ i1nu x 2mi y \ i1nu x 2mi y \\ 

y(u ,u } ) = y + yi cos — — + cos + — — + cos — 

1 v V3V v a V3V v a V3V' 

i . /4nu y \ . /2nu x 2mi y \ , ,2mi x 2nu y \\ 
+y 2 sin(^)-sm( + ^) + sinl — ) (3) 

The coefficients yo, yi and y2 can be determined by the fitting to the GSFE surface that can be obtained by 
our first-principles calculations, yo represents the ground state energy , which does not influence the shape of 



GSFE surface and can be chosen to zero, so only j\ and 72 need to be fixed. We also list the fitted parameters 
71 and 72 of Eq. (|3]> for all calculated FCC metals in Table 2. Eq. Q is accurate to describe translational 
and reflection symmetry of the GSFE surface on (1 1 1) plane in FCC metals, and facilitates the computation of 
dislocation properties. 

3 The dissociation of dislocation and the governing equation 

3.1. The improved Peierls-Nabarro dislocation equation 

In the (111) plane of FCC metals, the primary slip system is (112)(111), so dislocations with Burgers' 
vector ^(110) will energetically favorable dissociate into two Shockley partials connected by a stacking-fault 
(SF) ribbon according to Frank's rule [ 1 1. The dissociation process is described by the following reaction 

I[Tl ] -» I[Tl2] +SF + hill]. (4) 
2 o o 

The dissociation can be understood from the GSFE surfaces for slip plane (111). Fig.2 shows that the lowest 
energy path happens to be along the (112) direction, and perfect dislocations (Burgers' vector b = ^[110]) is 
likely to dissociate into partial dislocations (Burgers' vectors are bi = |[1 12] and b2 = ^[112]) due to the lower 
energy barrier. The two partial dislocations are 60° fractional dislocations viz. the included angle between the 
dislocation line and the Burgers vector is 60° . The magnitude of the total Burgers vector is b = \[aj2 and a is the 
lattice constant for FCC metal. The Burgers vectors of the two partial dislocations have magnitudes b\ = b2 = 
V3&/3. Hereafter, let the dislocation lines parallel to the y-axis ([1 12] direction) and perpendicular to the total 
Burgers vector of edge dislocation ([1 10] direction). The x-component and y-component respectively represent 
the edge component and screw component of partial dislocations for dissociation from edge dislocation. 

The 2D dislocation equation for straight dislocations that describe the balance of atoms on the border based 
on the lattice dynamics and the symmetry principle takes the following form lPT8l 

_p_ 1 <£u^_K 1 p° dx ,du x ^ = 
2 dx 2 2n J_ 0o x — x^ dx >\x=x' 

fodh^ K, f^_dx_^\\ pQftf) (5 ) 

2 dx 2 2n J_ 00 x — x^ dx "x=x 

where K e and K s are the energy factors of the edge and screw dislocations, u x and u y are the edge and screw 
components of displacements. For the isotropic solid, K e = /mt/(1-v) and K s = /ucr, with /u being shear modulus 
and v Poisson's ratio, o~ is the area of the primitive cell of the misfit plane. The shear modulus and the Poisson's 
ratio can be expressed as the second-order elastic constants c,y, i.e.,// = (c\\—c\<i)I2 andv = cnl(c\\ +C12) (the 
values as shown in Table 1). The coefficients of the second-order derivations fl e and fl s relate with the acoustic 
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phonon velocity and the lattice geometry structure 

3 2 2v 

Pe = 4^^^ - tan 2 cos 2 0), (6) 
3 

P s = -tan 2 6>sin 2 0), (7) 

where 9 and cf> are the orientation angles of the relative position of a pair of neighbor atoms in the intrinsic 
frame with the axes given by the polarization directions, and 9 = n/4 and <p = n/6 for FCC crystals, O is 
the volume of the primitive cell. Comparing with the generalized 2D P-N equation in P-N model, there are 
two extra second-order derivations that represent the discreteness effect of crystal in the new equation. From 
exactly solvable models iTTTl [TBI in which discreteness of a lattice can be fully considered, one clearly sees 
the discreteness correction appearing in terms of the second-differential of the displacement. Physically, the 
integro term in Eq. ([5]) represents a long range interaction which is inversely proportional to the distance, the 
differential term describes short range interaction which results from the interaction among the atoms on the 
surface plane. While f5 e and fi s are taken to be zero namely the discreteness effect is neglected the generalized 
2D P-N equation can be obtained [16]. 



3.2. The core properties of dissociated dislocation with variational method 

The dislocation solution of 2D dislocation equation is unknown. In this paper, the variational method 
is applied for the 2D dislocation equation lfT7l l35l . It can be straightforwardly verified that the variational 

functional of the 2D dislocation equation takes the following form 

-+°° Ai.x^ v r +°° r+°°,]u x du x x-r 

dxdx 



du x du x 


x=x 


x - X 


dx dx 


d 


du y du y 


> 

x=x 


X - X 


dx dx 


d 



-»+oo r*+oo 

/III- ill! \ — x 

dxdx 



" J —oo UA ^ n J — oo J— o 

X+oo 
y(u x , u y )dx, (8) 
oo 



where d is the distance of the nearest atoms in the misfit plane. 

In order to obtain the structure of dissociated dislocation, we represent each of the two partials by two 
closely spaced Peierls dislocations (of arctan-type) that can assure the correct asymptotic behavior. Thus, trial 
dislocation solution for dissociated dislocation can be written as 

b , ( x + d eq l2 ( x-d eq l2 M b 
u = — arctan ) + arctan + -, 

u y = ^jarctan( -^—)- arctan! j 3 —)), (9) 

where d eq is the equilibrium separation of two partial dislocations and ±d eq /2 are the positions of the two 
partials, g e and g s represent the core width of the edge and screw components respectively. The relation between 
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the edge and screw components of the relative displacement determines the dissociation path in the glide plane. 
Generally, the width of edge component is not equal to that of screw component and the dissociation path 
deviations from the crystallographic Burgers vector |[T6l . Thus, the dissociated dislocations cannot be dealt with 
ID dislocation equation. 

Substituting the trial solution into the variational functional Eq.([8]), the total variational functional / is as a 
function of d eq , £ e and £ s . The separation and the width of edge and screw components can be determined by 
solving the following equations 

dJ 

= 0, (10) 



dd eq 



£-a (id 

< 12) 

After a tedious although straightforward algebra one can obtain that the variational parameters of the trial 
dislocation solutions Eqs.Q for all our calculated FCC metals. Table 3 gives observed and calculated disso- 
ciated distances for the Al edge-character partial dislocation, in comparison with the wide range of previous 
results based on P-N model, atomistic simulations and density-functional theory (DFT). The partial separation 
distance in our results is 3.6ft (where ft is length of the total Brugers' vector), which is in good agreement with 
the experiment value 2.8ft obtained from weak-beam transmission electron microscopy [5] and the simulated 
value 2.4ft ~ 3.4b from DFT method in which the periodic supercell uses 552 atoms in a large radius cylindrical 
slab surrounded by vacuum H, as well as the value 3.2b obtained from the molecular dynamics simulation with 
the glue potential [10]. An earlier DFT calculation, which employed the quasi-continuum (QC) DFT method, 
where a small first principles cell containing the dislocation is embedded in the strain field produced by atoms 
interacting through an interatomic force law based on the embedded atom method (EAM), gives a smaller dis- 
sociated distance 2.0ft [7 ]. The difference of the QC-DFT-EAM calculation most likely stems from the small 
(84 atoms) supercell used for the first-principles calculations. The P-N model, which is a semidiscrete model, is 
essentially a continuum treatment. But the atomic-scale discreteness of crystals really matters in the region of 
partial dislocation core, so P-N models [7 , 8 ] are not well suited for describing details of the strain field in par- 
tial core. Different implementations of the P-N model predict different dissociated distances, though both are 
based on first-principles calculations of the GSFE surface. Several atomistic simulations for edge dislocation in 
Al have obtained different results for the dissociated distance between two partials, ranging from 1.7ft ~ 5.6ft 
ifTUl [TT1 [T2l : this variation is possible the result of different treatments for the boundary conditions 01. In our 
calculation, the GSFE surface is calculated from first-principle calculations accurately and the 2D improved 
P-N equation Eqs.(|5]> fully includes the discreteness effect of crystal, so our method is appropriate for dealing 



with the Shockley partials which are the mixture dislocations including edge component and screw component. 
The disregistry and the dislocation density which is defined as p = du/dx for Al, determined from our first 
principles calculations combined with the improved P-N equation, are shown in Fig. 3 (a) and (b), respectively. 
It is worth noting that the maximum value of relative displacements along y-axis u y for Al is 0.13ft, and is much 
smaller than the projection on the y-axis of the Burgers' vector of the partial g(H2) = 0.29ft (the black dot 
in Fig.3(a)). This means that the area within two partials is not a pure stacking fault ribbon, where the core 
regions of the two partials overlap and the displacement of y-component annihilate with each other due to the 
opposite direction 151 l37l. The structure of dislocations in Al cannot be dealt with the ID model due to the 
strong overlap of two partials as illustrated in Fig.3(b). The two partial dislocations can not be identified due to 
the large overlap that is consistent with the experiment and the first principle simulation |@). 

In this work, the core structures of edge dislocation in other FCC metals Cu, Ir, Pd, and Pt have also 
been calculated by using first principles calculations GSFE surfaces based on the improved P-N equation. The 
core widths both of the edge and screw components and the dissociated distances are listed in Table 4, and 
the disregistry and the dislocation density for Cu, Ir, Pd, and Pt are shown in Fig.4, Fig. 5, Fig. 6, and Fig.7, 
respectively. We also listed the results obtained from the P-N model in which the second-derivative term is 
neglected in Table 4. One can see that the discrete effects always increase the widths both for edge and screw 
components. Physically, the improved P-N equation Eq.Q is more functional to determine the structures of 
dislocation than the classical P-N equation [ 18 ], but for the calculated FCC metals we regret not being able to 
find any experimental data which would clearly demonstrate that the results with second-derivative term are 
more superior than those without second-derivative term. We cannot be sure that which model is a better one 
at the present level of knowledge. However, we believe the new experimental results will clarify this question 
mark. 

It has been known that the deformation mechanism cannot be explained by the absolute value of ISFE yisF 
alone and the entire behavior has to be understood in terms of the ratio JusfIjasf along (112) lfT4l . so we also 
list the ratios of 7usf/7isf in Table 4. When this ratio is low, then the energy barrier that has to be overcome for 
creating a trailing partial is very low and therefore it will be possible to observe the full dislocation. Pt has the 
lowest ratio 7usf/tisf value 1.10 in our calculation, so the overlap of two partials is so strong that one cannot 
distinguish them though dissociation behavior exists. In Fig .7, we can see that the maximum value of u y for Pt 
is 0.11ft and two-peak of the dislocation density for edge component p x vanishes. But when this ratio is large, 
the energy increase necessary for nucleating the trailing partial substantial, which is in the case of Cu with the 
large ratio 7usf/tisf value 4.07. In Fig. 4, we can indicate that the overlap of two partials is weak in Cu. It 
has verified in simulations (Ref[38|) extended partial dislocations in Cu have been observed as the predominant 
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deformation mechanism at nanocrystalline grains. It is worth noting that though the dissociated distance 5.5b 
for Pt is larger than that 5.4b for Pd, the core width of Pt is larger than that of Pd. This explains why Fig. 6 
for Pd shows two-peak of p x but Fig .7 for Pt does not. Our detailed calculations demonstrate that whether full 
dislocations or Shockley partials will easily be observed in FCC lattice must be understood in terms of the ratio 
Tusf/Tisf- 

4 conclusion 

In conclusion, we have performed first-principles calculations to obtain the GSFE surface for the { 1 1 1 } glide 
plane in FCC metals Cu, Al, Ir, Pd, and Pt. The accuracy of the method have been tested by calculating values 
for various stacking fault energies which favorably compare with the previous theoretical and experimental 
results. From these calculations, we extract the core properties for the edge dislocation of dissociation into 
partials, using the improved P-N equation which includes the discreteness effect of crystal. The Ritz variational 
method is presented to solve the dislocation equation and the trial solution is constituted by two arctan-type 
functions which represents the dislocation dissociating into two partials. The core structures, including the core 
widths both of the edge and screw components, dissociation behavior for edge dislocations have been inves- 
tigated. The dissociated distance between two partials for Al in our calculation agrees well with the values 
obtained from the numeric simulation with DFT and MD methods, as well as experiment. Our detailed calcu- 
lations demonstrate that whether full dislocations or Shockley partials will easily be observed in FCC lattice 
must be understood in terms of the ratio 7usf/7isf, and show that it is preferred to create partial dislocation in 
Cu, and to be easily observed full dislocation in Al, Ir, Pd, and especially Pt 
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Table 1 The calculated equilibrium lattice parameters a [A] and elastic constants [GPa], comparison with the experimental and previous 
calculated results. 





a 


Cll 


Cl2 


c 44 




V 


Al 


4.05, 4.05" 


112.3, 114.3* 


60.7,61.3'' 


30.2,31.6* 


25.8, 26.3* 


0.350, 0.349* 


Cu 


3.64, 3.62" 


175.3, 176.2* 


124.4, 124.9* 


80.9, 81.8* 


25.5,25.7* 


0.415, 0.415* 


Ir 


3.88, 3.84° 


579.6, 582.3* 


240.8, 241.3* 


261.2, 262.0* 


167.9, 170.5* 


0.294, 0.293* 


Pd 


3.96, 3.89" 


232.1, 234.1* 


175.2, 176.1* 


70.5, 71.2* 


28.5, 29.0* 


0.430, 0.420* 


Pt 


3.99, 3.92" 


356.8, 358.0* 


252.2, 253.6* 


76.3, 77.4* 


52.3, 52.2* 


0.414, 0.415* 



"Referenced 



*Referencet30l 



Table 2 The various stacking-fault energies on (1 1 1) plane for FCC metals. The experimental and the other calculated results are also 
listed. All values are in units of ml/m 2 . 





Cu 


Al 


Ir 


Pd 


Pt 




-160 


-158 


-615 


-167 


-171 


72 


-260 


-213 


-928 


-243 


-188 


tisf 


43 


158 


359 


122 


282 


yisF 11 


75 


146 


41 


244 


356 


risF 6 


43 


150 




101 


111 


7isf c 


51 


153 




186 




7isF rf 


40, 45, 169 


166 


300 


180 


322 


yusF along <112> 


175 


225 


753 


215 


311 


7usf along (110) 


638 


633 


2462 


668 


685 



"Reference|31 1, by EAM calculations 
*Reference| 32|, by MEAM calculations; 
c Reference ll3l . by first-principles calculations; 
rf Reference[ 1 33 34], by experiments. 



Table 3 Shockley partial dissociated distance d eq for 1/2(110) edge dislocations in Al. Our calculated result is compared with the 
experimental results , the previous theoretical values including the Peierls-Nabarro methods and the numerical simulations. All data 
are in units of the length of Burgers vector b = y/2/2a. 



This work a Ref|5| 


b Ref[4] b Ref|6| 


c Ref|7| 


c Ref| 8 1 


d Ref(9) 


d Ref|10| d Ref|ll| 


3.6 2.8 


2.5 ~ 3.4 2.0 


1.2 


2.6 


3.2 


5.1 5.6,1.7-5.2 



" experiment from weak-beam transmission electron microscopy (WB-TEM) 
* obtained from density-functional theory (DFT); 
obtained from Peierls-Nabarro model; 
'' obtained from atomistic simulations 



Table 4 The predicted core structures, including core widths both the edge and screw components, dissociated distance in FCC metals 
Cu, Al, Ir, Pd, and Pt are listed. For completeness, we also list the results obtained from the P-N model in which the second-derivative 
term is neglected. For discussion, we give the ratio 7usf/7isf values along the (112) direction. The ratio of 7usf/Tisf is dimensionless, 
and £ s and d eq are in units of the length of Burgers vector b = V2/2a. 





Cu 


Al 


Ir 


Pd 


Pt 




2.6°, 2.1* 


2.5°, 2.0* 


3.5", 3.1* 


3.1", 2.4* 


4.8", 4.0* 


e, 


1.9°, 1.5* 


2.1°, 1.6* 


2.9", 2.5* 


2.4", 1.8* 


3.9", 3.0* 




5.5", 3.7* 


3.6", 2.6* 


5.6", 4.6* 


5.4", 3.5* 


5.5", 4.0* 


Tusf/Tisf 


4.07 


1.42 


2.10 


1.76 


1.10 



" obtained from improved P-N equation (with second-derivative term) 
* obtained from P-N equation (without second-derivative term). 
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Fig. 1 The GSFE surface for displacements along a (111) plane in Al (the corners of the plane and its center correspond to the 
configurations of ideal Al lattice). The first energy maximum encountered along the (112) direction is the USFE; the first energy 
minimum value at 1/6(112) corresponds to the ISF configuration. The projection of the GSFE surface along the 1/2(110) is symmetric 
with respect to the slip displacement of V2a/4. 



800 




6 0.2 0.4 0.6 0.8 1 

Slide along 1/6<112> 



Fig. 2 Projection of the y-surface in the [112] directions on the (111) plane. There have been termed the intrinsic stacking fault (ISF) 
energies corresponding to the final positions and the unstable stacking fault (USF) energies corresponding to the maximum of the 
curves. 
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Fig. 3. Disregistry profile u(x) (a) and dislocation density (b) for the edge dislocation in Al are obtained from our first-principles 
calculations combined with the improved P-N equation. The solid line and the dashed line represent the edge and screw components 
of displacements, respectively. The black dot marks the displacement of the screw component of the crystallographic Shockley partial 
dislocation, which is about 0.29b. Here the maximum value of relative displacements of screw components of the partial dislocation is 
0.13b. 





Fig. 4. Same as in Fig. 3, but for Cu. Here the maximum value of relative displacements of screw components of the partial dislocation is 0.20b. 
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Fig. 7. Same as in Fig. 3, but for Pt. Here the maximum value of relative displacements of screw components of the partial dislocation is 0.1 lb. 
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